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ABSTRACT 

We present a period-luminosity-amplitude analysis of 5899 red giant and binary stars in the Large 
Magellanic Cloud, using publicly available observations of the MACHO project. For each star, we 
determined new periods, which were double-checked in order to exclude aliases and false periods. The 
period-luminosity relations confirm the existence of a short-period, small-amplitude P-L sequence at 
periods shortward of Seq. A. We point out that the widely accepted sequence of eclipsing binaries 
between Seqs. C and D, known as Seq. E, does not exist. The correct position for Seq. E is at periods 
a factor of two greater, and the few stars genuinely lying between Seq. C and D are under- luminous 
Mira variables, presumably enshrouded in dust. The true Seq. E overlaps with the sequence of Long 
Secondary Periods (Seq. D) and their P-L relation is well described by a simple model assuming 
Roche geometry. The amplitudes of LSPs have properties that are different from both the pulsations 
and the ellipsoidal variations, but they are more similar to the former than the latter, arguing for 
pulsation rather than binarity as the origin of the LSP phenomenon. 

Subject headings: binaries: eclipsing — galaxies: individual (Large Magellanic Cloud) — stars: AGB 
and post-AGB — stars: oscillations — stars: statistics — stars: variables: other 
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1. INTRODUCTION 

The multiplicity of red giant period-luminosity (P-L) 
relations has been a major discovery on the road to inter- 
preting complex light vari ations of these stars . Fol lowing 
the tw o seminal papers bv lWood et alJ l|1999l ) and I Wood! 
( 2000), a picture has emerged that can be summarized 
as follows: large-amplitude Mira stars pulsate in the fun- 
damental mode, whereas smaller-amplitude semiregulars 
are often multimode pulsators, in which various overton e 
modes can be excited (see also lBedding fc Ziilstraf 19981 . 
Besides th e pulsating P-L seq uences (Seq. A, B and C, as 
labeled bv lWood et alJ l)1999t) 1. two other sequences were 
suggested: Seq. E with red giants in eclipsing binaries 
and Seq. D with stars that have long secondary periods 
(LSPs). The latter pose a great mystery and the na- 
ture of their slow variations is still not understood, with 
several different mechan isms proposed l|01ivier fc WoodI 
20031 'Wood et a,Ul200l . 

The basic picture of multiple P-L relations has been 
confirmed by many independent studies, mostly based 
on if-band magnitudes. It has emerged that the origi- 
nal five sequences have further details, including a break 
at the tip of the Red Giant Branch (RGB), which is 
due to the existence of distinct RGB pulsators that 
are mixed with the more evolved AGB vari a bles (e.g. 
Ita et alJ l200l IKiss fc Bedding] 12001 12001 lit a et all 



2004albUSoszvnski et al.ll2004aUFraser ■ et al]f2005|) . Al- 



most all authors have accepted the existence of the dis- 
tinct sequences of red giant binaries (S eq. E) and LSP 
stars (Seq. D). The only exception was iSoszvhski et alJ 
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( 2004b]), who showed that Seqs. E and D seem to merge 
at a specific luminosity (as measured by the Wesenheit 
index) and suggested that this may imply the binary ori- 
gin of LSPs. 

Here we report on a combined analysis of MACHO 
observations of eclipsing binaries and red giants in the 
Large Magellanic Cloud, which shed new light on these 
stars and on the LSP phenomenon. 

2. DATA ANALYSIS AND RESULTS 

Our results are based on two sets of publicly avail- 
able almost eight-year long MACHO light curves. 
A detailed d escript i on of the MACHO project can 
be found in ICookl lfl99^ . Some of the data are 
offered for download through the MACHO website 
(http:/ /wwwmacho.mcmaster.ca), where one can choose 
specific samples based on an automated classification of 
variability type. Using the Web interface, we individ- 
ually downloaded all light curves classified as eclipsing 
binaries (6833 stars) and as red giant variables (classi- 
fied as Wood A, B, C and D classes; 2868 stars). 

In the case of the eclipsing binaries, it became obvi- 
ous very quickly that the classification was not perfect, 
and a large fraction of stars turned out to be Cepheids, 
RR Lyrae stars or long-period variables. We also found 
that the catalogued periods were incorrect for a signifi- 
cant number of stars. We therefore reclassified all 6833 
stars and re-determined their periods, using the following 
procedure (more details will be given elsewhere). 

Periods were first estimate d using the Phas e Dis- 
persion Minimization method ( Stel lingwerJ n.978) . We 
then checked all the folded light curves by eye and 
refined the periods w i th the Strin g-Length method 
l|Lafler fc Kinmanlll965l lClarkdl2002T) . which is more re- 
liable than PDM when the light curve contains long flat 
sections and very narrow minima, as is the case for many 
eclipsing binaries. Also, in many cases PDM gave har- 
monics or subharmonics of the true period, which was 
only recognized through the visual inspection of every 
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phase diagram. We also examined the color variations 
to identify and exclude pulsating stars with sinusoidal 
light curves. After this analysis, 3031 stars remained as 
genuine eclipsing or ellipsoidal variables. 

Next, we classified the binary sample using Fou r ier de- 
composition of their phase diagrams. Rucinski (1993) 
showed that light curves of W UMa systems (contact bi- 
naries) can be quantitatively described using only two 
coefficients, &2 and en, of the cosine decomposition 
a-i cos(27ri</?). iPoimahskl l)2002[) tested the behavior of 
semi-detached and detached systems in the &2 — &4 plane 
by decomposing theoretical light curves into Fourier co- 
efficients. We found that only stars with "W UMa-like" 
light curve shape composed the sequence (which is plot- 
ted in Fig. ^ , while detached and semi-detached systems 
are spread everywhere in the P-L plane. 

Our second set of light curves were those of the 2868 
publicly available MACHO red giant variables. Since 
they often show multiply periodic light variations, we 
determined periods with iterative sine wave fitting. As a 
measure of significance, we also estimated the S/N ratio 
of the peaks in the Fourier spectra. Since the noise in 
the Fourier spectrum increased toward lower frequencies, 
different values of the S/N were used when determining 
whether a peak was real for different period-luminosity 
relations (S/N cutoff was set to 3 for Seq. A', 4 for Seq. 
A, 6 for Seqs. B, C and 10 for Seq. D). We omitted 
periods close to 1 yr, because many light curves show 
variations with this period that are not real. As a result, 
a total of 4315 significant frequencies were identified for 
the 2868 stars. 

We studied these two samples in the P-L plane. In 
order to reduce the effects of interstellar extinction and 
allow a direct comparison with previous results in the 
literature, we plotted the period-fC magnitude rela- 
tion. We obtained near infrared magnitudes by cross- 
correlation with the 2MASS All-Sky Point Source Cata- 
log (http://irsa.ipac.caltech.edu), with a search radius of 
3". The resulting P-L diagram is shown in the left panel 
of Fig. [U 

To our surprise, the P-L relation of the binary sample 
did not follow Seq. E, as we had expected. Instead, they 
overlapped with Seq. D, which at first sight appears to 
give strong evidence for the binary origin of Seq. D and 
prompted us to investigate the issue in more detail. 

3. DISCUSSION 

The combined P-L plot in Fig. ^ shows the well- 
known complex structure of distinct sequences. Besides 
sequences A, B, C, D and E, we also detect the ex- 
istence of the faintly visible new short-period P-L se- 
quence (Seq. A') on the left-hand si de boundary of 
the diagram (label ed as m an d b4 i n iSoszvhski et alJ 
l|2004a|) and P 4 in IKiss fc Labi p006j) L The eclipsing 
stars (pluses) seem to merge with Seq. D rather than 
forming Seq . E as adopted in the literature llWood et al. 
[19991 | Woodl200?i IKiss fc Beddindl2003t 12004 llta'et'aL 
l2004alrtlNoda et al .ll2004UFra,ser et al .11200,51) . 

To clarify this issue, we r e-checked per iods for: (i) stars 
on Seq. E in Fig. 1 of iWoodl (|2000|) . for which the 
identifiers and basic data were kindly provided by Peter 
Wood; (ii) stars on Seq. E in our Fig. (black dots). 
In both cases, it turned out that for most of the objects, 
the given periods were half of the true ones, as one might 



expect from a Fourier analysis of eclipsing binary light 
curves. We have carefully double-checked all individual 
light curves on Seq. E and D and corrected the periods 
(see Fig. 0). The final P-L plot is shown in the right 
panel of Fig. ^ 

As a result of the period correction, the sequence be- 
tween C and D, which is known as Seq. E in the liter- 
ature, has completely disappeared. A few stars remain 
in the gap but practically none of them are eclipsing bi- 
naries. The majority turn out to be Mira stars that are 
very red (J — K > 2 mag) and are presumably carbon- 
rich Miras that are dimmed by circumstellar dust clouds 
(bottom panel of Fig. EI) • W° propose to retain the label 
E for the sequence of ellipsoidal variables in its corrected 
position (right panel of Fig. U) . 

In the literature, only ISoszvhski et alJ l)2004bl) plotted 
the eclipsing binary sequence at the correct (doubled) 
period. However, they did not discuss this issue and one 
can still find more re cent studies wher e Seq. E was shown 
at the wrong period. ISoszvhski et alJ l)2004b() have, how- 
ever, shown that Roche geometry gives a good fit to the 
OGLE ellipsoidal variables. We have also checked this 
on the MACHO sample. We calculated the theoretical 
orbital periods of systems at mass ratios of 1, where the 
components fill their Roche lobes. For this, we used evo- 
lutionary models of lCastel lanLe±_all 12QQ3 D and applied 
equations of Section 4 in lSoszvriski et alJ l)2004bfl . (Note, 
however, that the definition of f(q) in their Eq. 1 actu- 
ally gave f(q)~ 1 as the filling factor.) For the calculations 
we took the evolutionary tracks of the 0.85 M Q and 2.5 
Mq models, since theses masses represent the mass limits 
of stars that evolve through the RGB and AGB. The K 
magnitudes of the models were deter mined from T P ff vs. 
(V — K) c alibrations, comb ined from iHoudashelt et alJ 
•21)00) and lKucinskas et all l)2006[) . 

In the right panel of Fig. ^ we show these two lim- 
its (the shorter period line belongs to the higher mass). 
For any other mass ratios, the orbital periods shift to- 
wards smaller values. This simple approach with Roche 
geometry describes remarkably well the observed period- 
luminosity relation of ellipsoidal variables in the MACHO 
sample , in agreement with the study of ISoszvhski et aD 
(2004b). However, the OGLE sample contains a larger 
fraction of ellipsoidals and most of them have longer pe- 
riods than are predicted by the models. Those stars d o 
not entirely fill the Roche lobe ( Sos zvhski et al.l l2004b) . 
It is also worth mentioning that the low fraction of cllip- 
soidal /eclips ing RGB sta r s in t he MACHO data (1.5%) 
was used by iWood et alJ l)2004|) as an argument against 
the binary origin of LSPs. However, the OGLE statistics 
clearly showed that there are at least 10 times more such 
RGB stars, and so that argument is no longer valid. 

Since the lower mass model fits Seq. D quite well, 
the question arises: how similar are the ellipsoidal and 
LSP variables? To assess this, we examined the ampli- 
tudes of these stars. ISoszvhski et alJ l)2004b|) mentioned 
that amplitudes of LSPs are positively correlated with 
the brightness of the star. Compared to OGLE data, 
the MACHO observations have the advantage of giving 
information on the color variations. We examined the 
amplitudes in the MACHO blue and red bands of the 
binary and LSP stars, and also included Seq. C stars, 
which allow a comparison with stars that we know to be 
pulsating. 
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Fig. 1. — Period-Luminosity relations for red giant pulsators (black dots) before (left panel) and after (right panel) period correction of 
stars on Seq. E; 310 ellipsoidal and contact eclipsing binaries (red pluses) are shown with the correct periods in both plots. The solid lines 
in the right panel show model calculations using evolutionary tracks and assuming Roche-geometry (see text). 
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Example of phase diagrams before (left panels) and 
after (right panel) period correction. The first and the second 
star from the top show eclipsing binary and ellipsoidal variable 
light curve, while the third star has a typical LSP variation. The 
bottom panel shows the light curve of a Mira star that lies in the 
gap between Seqs. C and D and which faded ~2 mags during 8 
years of MACHO observations. 

The amplitudes were measured by fitting smooth spline 
functions to the phased light curves. The resulting blue 
peak-to-peak amplitudes are shown as a function of K 
magnitude in Fig. [3 where several features are apparent. 
Firstly, the amplitudes of ellipsoidal variables (circles in 
the upper panel) do not show any correlation with lumi- 
nosity, as expected for variation caused by the geometry 



Fig. 3. — Luminosity-amplitude relation of stars on Seq. C, 
D and E. Top panel: Black crosses show the luminosity-amplitude 
relation for Seq. D stars, while turquoise circles the same for Seq. 
E stars, which contain ellipsoidal variables. The difference of the 
two sequence is very clear, e.g. amplitudes of Seq. D stars show 
moderately positive correlation with the luminosity. Bottom panel: 
Luminosity-amplitude relation for Seq. C stars, which are funda- 
mental mode pulsators. This relation is significantly different from 
those of Seqs. D and E stars. 

of a binary system. Secondly, the amplitudes of LSP vari- 
ables (crosses in the upper panel) increase with luminos- 
ity, with their distribution forming a striking triangular 
envelope (the correlation coefficient is r ~ 0.51). It is also 
very interesting that the points below the well-defined 
upper envelope seem to show a flat distribution. What- 
ever the cause of the LSP phenomenon, there appears 
to be a maximum possible amplitude at each luminosity, 
with an apparently uniform distribution of amplitudes 
beneath this maximum value. 

In comparison, stars on the pulsating P-L sequence C 
(lower panel in Fig. |3J), show a very different distribu- 
tion. To quantify the difference between the amplitude- 



4 



Derekas et al. 




log P (d) 



Fig. 4. — Blue-to-red amplitude ratios as function of period 
for stars on Seqs. B (grey triangles), C (orange pluses), D (black 
crosses) and E (purple circles). Grey triangles show those short 
period pulsations of Seq. D stars, which fall to Seq. B. The two 
clearly pulsating sequences, Seqs. B and C are the continuation 
of each other. Ellipsoidal variables (Seq. E) have lower amplitude 
ratio than pulsating stars. Amplitude ratios of Seq. D stars are 
more similar to those of pulsating stars. 

luminosity distribution of S eqs. C and D, w e perfo rmed 
the multivariate e-test of iSzekelv fc Rizzol l)2004D . ap- 
plied for two dimensions. In this test we measure the 
effects of random permutations of the initial distribu- 
tions via changes in the "information energy" of the two 
distributions. In our case, 1000 random permutations 
showed that the difference between the C and D samples 
is highly significant. Therefore, we conclude that the 
physical mechanism causing the LSP phenomenon must 
be different from both the ellipsoidal variations of Seq. E 
and the radial fundamental- mode pulsations of Seq. C. 
Note, however, that for Seqs. A and B, there is a similar 
positive correlation between amplitude and luminosity 
(see th e up per three panels i n figs . 4 in lKiss fc Be dding 
( 2003| and lKiss & Beddind pOOl l than for Seq. D. 

At the same time, comparing blue and red amplitudes 
for individual stars revealed further interesting informa- 
tion. In Fig. 0] we plot the blue-to-red amplitude ratios 
as function of period for Seqs. B, C, D and E stars. For 
the pulsating objects the median of the ratio is 1.40, in- 
dicating strong color, thus temperature, changes during 
the pulsations. For ellipsoidal/eclipsing binaries, the me- 
dian ratio is 1.13, while the LSPs have a median ratio of 
1.29, being more s imilar to the pulsating stars (see also 
iHuber et al1l2003D. This behavior agrees with the find- 
ings of lWood et all l)2004fl who, based on color-amplitude 
variations of single objects, argued for the pulsational ori- 
gin of LSPs. The overall statistics of more than 700 LSP 
variables favors this arg ument over t h e bina ry hypoth- 
esis, also agreeing with IHinkle et alJ lj2002f) . who con- 
cluded that the long-period velocity changes in their ob- 
served stars probably result from some kind of pulsation. 
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4. SUMMARY 

The main results of this paper can be summarized as 
follows: 

• the period-luminosity relations of ^6000 stars 
based on MACHO data confirm the existence of 
the short-period, small-amplitude P-L sequence at 
shortward of Seq. A, which belongs to a higher- 
overtone pulsation mode. We label this Seq. A'. 

• the widely accepted sequence of eclipsing binaries 
between C and D, known as Seq. E, does not exist. 
The correct position for Seq. E, which comprises 
contact binaries and ellipsoidal variables, is at pe- 
riods a factor of two greater. The true Seq. E 
overlaps with the LSPs (Seq. D), which appears to 
suggest a binary origin for the LSP phenomenon 
(but see the last point). 

• of the few stars that genuinely lie between Seq. C 
and D, most are under-luminous Mira variables, 
presumably enshrouded in dust. 

• we confirmed that ellipsoidal variables have a simi- 
lar P-L relation to LSP stars. Their P-L relation is 
well described by a simple model assuming Roche 
geometry. 

• the amplitudes of LSPs have properties that are dif- 
ferent from both the pulsations and the ellipsoidal 
variations, but they are more similar to the former 
than the latter, arguing for pulsation rather than 
binarity as the origin of the LSP phenomenon. 
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